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SUMMARY 


Results of tests of a 35^ swept-wing fighter airplane, during which 
lateral oscillations were performed over a Mach number range from 0.^1 to 
0.79 at a pressure altitude of 10,000 feet and from 0.49 to 1.04 at 35^000 
feet, are presented in this report. Experimental and computed values for 
the period of the lateral oscillation and time required to damp to half 
amplitude are shown. One sample oscillation time history is included for 
each test altitude. 

The airplane was found to be laterally stable, statically and dynami- 
cally, throughout the range of speeds tested. At both altitudes, the 
variation with Mach momber of the period of the lateral oscillation was 
satisfactorily predicted from available and estimated aerodynamic and mass 
parameters. The time required to damp to half amplitude, as measured from 
flight at both altitudes, varied with Mach number in essentially the same 
manner as predicted from computations. The measured damping was somewhat 
better than that obtained from computations for the altitude of 35^000 
feet, particularly at a Mach number of O. 92 . An increase in time to damp 
to half amplitude was noted between Mach numbers of 0.95 Q-^d l.OU. 


INTRODUCTION 


As part of a general research program concerned with the lateral 
dynamic stability and handling characteristics of high-speed, high-altitude 
airplanes, the Ames Aeronautical Laboratory of the NACA has tested a 35 ^ 
swept-wing fighter airplane through a wide range of flight speeds and 
altitudes . 

The purpose of this report is to present results of tests of the 
lateral oscillatory characteristics made during a series of four flights. 
Comparisons are included of the computed variation of period and damping 
of the lateral oscillation with the measured values. These comparisons 


^Supersedes recently declassified NACA RM A 5 IC 28 by Walter E. McNeill 
and George E. Cooper. 
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indicate the accuracy with which the oscillatory behavior of an airplane 
can be predicted under various flight conditions using available or esti- 
mated mass parameters and stability derivatives, and neglecting such 
effects as aeroelasticity and iinsteady lift. 


SYMBOLS 


Cl 

Cy 

Cy 

CYt, 
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lift coefficient, 
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lateral-force coefficient, 
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— per radian 
^2V 

moment of inertia about flight-path axis, slug-ft^ 

moment of inertia about axis normal to flight path in the plane 
of symmetry, slug-ft^ 

Mach nimber 

period of oscillation, sec 
wing area, sq ft 

time to damp to half amplitude, sec 

true airspeed, ft/sec 

weight of airplane, lb 

wing span, ft 

pressure altitude, ft 

distance parallel to longitudinal reference axis from center of 
gravity of the airplane to center of pressirre of vertical tail 
in yaw, ft 

rolling angular velocity, radians/sec 

dynamic pressure, Ib/sq ft 

yawing angular velocity, radians/sec 

normal distance from longitudinal reference axis to center of 
pressure of vertical tail in yaw, ft 

angle of attack of longitudinal reference axis, deg 

angle of attack of longitudinal reference axis for zero lift, deg 

angle of sideslip, radians 

dihedral angle, radians 

angle between longitudinal reference axis and principal axis of 
airplane, positive when reference axis is above principal axis 
at nose, deg 
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a sidewash angle at vertical tail resulting from the wing in rolling 

flow^ positive for positive lateral force, radians 


tp 

|P| 


angle of bank, radians 

ratio of amplitude of angle of bank to amplitude of sideslip angle 
for the oscillatory mode 


Subscripts 


h contributed by horizontal tail 

i+f contributed by fuselage and wing-fuselage interference 
M pertaining to a given Mach number 

V contributed by vertical tall 

w contributed by wing 

W.T. obtained from wind-tunnel tests 


INSTRUMENTATION AND FLIGHT TECHNIQUE 


The general arrangement of the test airplane is shown by a photograph 
(fig. 1) and a two-view drawing (fig. 2). The principal dimensions are 
listed in table I. 

Standard NACA recording Instruments were used to measure angle of 
sideslip, rolling and yawing velocities, pressure altitude, and airspeed. 
Aileron and rudder deflections were recorded by NACA instruments as well 
as on separate channels of a 36-channel oscillograph. The rudder deflec- 
tion was known to an acc\aracy of 0.1°, while the aileron deflection was 
known within 0.3°. The nose-boom airspeed system described in reference 1 
was used to detennine Mach ntunber and the static and dynamic pressure. The 
records were synchronized by a l/lO-second instrument timer . 

At a pressure altitude of 10,000 feet, lateral oscillation maneuvers 
were performed through a range of Mach numbers from 0.4l to 0.79* At an 
altitude of 35,000 feet, oscillations were performed at Mach numbers from 

0.49 to 1.04. 

All oscillations performed at 10,000 feet were excited by retiirns 
from steady sideslips. At 35^900 feet, the airplane was disturbed both 
by returns from steady sideslips and by abruj)tly deflecting the rudder 
and returning it to neutral, except at Mach numbers above 1.02, where 
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rudder kicks alone were used. During all test runs below a Mach number 
of 1.0, the rudder and ailerons were held essentially fixed following 
their return to trim positions with the aid of chains which prevented the 
pilot's moving his controls beyond a predetermined point. At Mach numbers 
above 1.0, chains were used on the rudder pedals only. 

All oscillations were performed in the clean condition and in level 
flight, with the exception of those at Mach numbers above 0.92 where dive 
angles up to 36° were required. 


RESULTS AND DISCUSSION 


Typical time histories of indicated airspeed, press\ore altitude, side- 
slip angle, rolling velocity, yawing velocity, total aileron deflection, 
and rudder deflection are shown in figure 3 for a Mach mmiber of O.79 at 
a pressure altitude of 10,000 feet. Figure 4 presents time histories of 
the same quantities for an average Mach number of 1.04 at about 35^000 
feet . 


The results of data obtained during similar lateral oscillations at 
altitudes of 10,000 and 35,000 feet are summarized in figure 5 in the form 
of period and time required to damp to half amplitude expressed as fimc- 
tions of Mach number. 


For comparison with the experimental results, ciirves of computed 
values for period and damping also are shown in figure 5* These values 
are solutions to stability quart ics derived from the lateral equations of 
motion presented in reference 2. The mass distribution and dimensional 
data used in computing period and damping were furnished by the manufac- 
turer. The methods used to measure or estimate the variation of the sta- 
bility derivatives with Mach number and lift coefficient are summarized 
in the appendix. All lateral derivatives were corrected for compressi- 
bility effects according to the Prandtl-Glauert rule, as outlined in the 
appendix, from M = 0 to M = 0.9. Each derivative was then plotted as a 
function of Mach nvimber-and the resulting curve was extrapolated at a con- 
stant slope from M = 0.9 to M = 1.0. Table II presents the values of the 
parameters used in computing period and damping at each Mach nvunber con- 
sidered at altitudes of 10,000 and 35,000 feet, together with the resulting 


values of P, T^/a^ and 


ItI 

|P|' 


The lift coefficients shown in table II are 


representative of the flight values except for Mach numbers greater than 
0.95^ where a deviation of less than I5 percent would be expected. Because 
of the small range of angle of attack involved (less than 5°)^ the rolling 
and yawing moments of inertia were assumed constant at the values given 
beneath table II. 


Figure 5 indicates that reasonably close correspondence (within 8 
percent) was obtained between the variation with Mach number of computed 
and measured values of period at pressure altitudes of 10,000 and 35^000 
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feet. The measured period is observed to decrease less rapidly with 
increasing Mach ntunber than did the computed value at 35 ,000 feet . At 
10,000 feet the opposite trend is seen; that is, the experimental value 
of the period decreased slightly more rapidly than did the computed value 
as Mach number was increased. No explanation for this phenomenon is 
apparent . 

Close agreement (within 7 percent) existed between the measured and 
computed values for T^/2 3-t Mach numbers below 0.6 at 10,000 feet. Above 

M = 0.6 the flight-test values for T1/2 became Increasingly higher than 
the computed values as the Mach number was increased. At 35,000 feet the 
meas\ired values for Tx/2 were somewhat lower than the computed values 
throughout the major portion of the Mach number range tested (0.U9 to 
1.04), with the best agreement occurring at low speed. Above M = 0.8 the 
experimental value of T1/2 decreased more rapidly, reaching a maximum 
deviation of about 20 percent from the computed curve at M = 0.92, then 
changed its slope gradually from negative to positive up to M = 1.04. 

Due to the scatter of test points at Mach nvunbers above 1.0, it is diffi- 
cult to determine more than the general trend indicated in figure 5* 

It is evident that a good prediction was made of the lateral period 
and damping of the test airplane for the range of lift coefficients con- 
sidered in the computations (O.O76 to 0.412) using Cj^p and C^p for the 
wing alone, as shown by the dashed curves in figure 5* this instance, 
nothing was gained by considering the contributions of the vertical tail 
in addition to the wing, as shown by the low and high Mach number points 
for which C^p and were computed by the methods of reference 3« As 

the lift coefficient of the airplane is increased, however, the vertical- 
tail contribution to Cnp becomes quite large and could be Included as 
shown in the appendix. 

Figure 6 presents the above flight information as the relationship 
between period and time to damp to half amplitude for each altitude , 
together with the corresponding computed values. Good agreement between 
the measiired and computed period-damping relationships is again demon- 
strated by this method of presentation, particularly for the altitude of 
10,000 feet below a Mach number of 0.6. As in figure 5, figure 6 shows a 
lower experimental value of T1/2 for 35,000 feet at all Mach numbers, 
particularly at M = 0.9, in relation to the computed values. 

The Armed Services lateral-oscillation specification which applied at 
the time of these tests (ref. 4 ) is shown in figure 6 for comparison with 
the characteristics of the test airplane. It is shown that the lateral 
oscillatory characteristics of the test airplane at 35,000 feet were 
entirely within the ^satisfactory region defined in reference 4 . The 
same figure shows that the airplane, at 10,000 feet, exhibited borderline 
characteristics with respect to the requirements of reference 4 except at 
Mach numbers between O.54 and O.79 where the period -damping relationships 
were in the satisfactory region. The computations indicate that for low 
Mach numbers, near O.35, characteristics exist which are satisfactory under 
the requirements of reference 4 . 
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According to the pilot’s comments, the lateral oscillatory charac- 
teristics of the test airplane were generally satisfactory at 10,000 feet. 
At 35,000 feet, the oscillations were somewhat ohjectionahle, partly 
because of the increased rolling that was present (see table III) and 
partly because of the noticeably decreased damping which was especially 
apparent in rough air. Rough air tended to prolong the oscillations at 
10,000 feet as well but, since there was considerably less rolling present 
at comparable Mach numbers, the motions were not considered so objection- 
able as those experienced at the higher altitude. 


Results of other lateral flying qualities investigations conducted 
at the Ames Aeronautical Laboratory (see ref. 5) have Indicated the pos- 
sibility of the use of the ratio of amplitude of angle of bank to amplitude 

I I 

of sideslip angle for the oscillatory mode as a criterion for satis- 

factory lateral oscillatory characteristics of fighter-type airplanes. 


Measured values of are presented, in addition to period and 

time to damp to half amplitude, in table III for the Mach niimber ranges 
covered at the test altitudes of 10,000 and ^'^,000 feet. 


Computed values of 


|cpl 

31 


are shown in table II for Mach numbers con- 


Through the speed ranges covered at 
both altitudes, for the test airplane was predicted with less accu- 


sidered at 10,000 and 35^000 feet, 
both altitudes, for the test 
racy than were period and damping. 


CONCLUSIONS 


1. In general, the lateral oscillatory characteristics of the test 
airplane were closely predicted from information based on wind-tunnel 
tests, although unsteady lift and aeroelastic effects were neglected. 

2. Throughout the range of Mach numbers tested (O.Ul to l.OL) the 
airplane was laterally stable both statically and dynamically. 

3 . The period of the lateral oscillation varied smoothly with Mach 
number over the range tested and was adequately represented by computed 
values at both test altitudes, with no error greater than 8 percent. 

L. The time required for the lateral oscillation to damp to half 
amplitude at test altitudes of 10,000 and 35,000 feet decreased with Mach 
number in essentially the same manner as indicated by computations, except 
at 2^,000 feet where the measured value of T 1/2 began to increase with 
Mach number above M = 0.95 • 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 28, 1951 
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APPENDIX 

ESTIMATION OF LATERAL- STABILITY DERIVATIVES 
FOR THE TEST AIRPLANE 


The values presented in table II for the lift-curve slope, angle of 
attack at zero lift, and the static lateral- stability derivatives were 
obtained from wind-tvinnel or other force-test methods and corrected for 
compressibility effects where tests did not cover the Mach number range 
considered in this report. The rotary derivatives were estimated by pub- 
lished theoretical methods applicable to swept wings. 

The procedures used in determination of the aerodynamic parameters 
and stability derivatives considered in this report are presented below. 


Longitudinal Trim Parameters 


Variation of lift-curve slope, with Mach number was determined 

from the results of subsonic tests in the Ames l6-foot wind tunnel and the 
Southern California Cooperative Wind Tunnel to a Mach number of 0.94 and 
supplemented by transonic btimp tests to a Mach number of I.06. 

Angle of attack for zero lift, ao^ was taken from unpublished results 
of NACA wing- flow tests. 


Static Lateral -Stability Derivatives 


Lateral force due to sideslip .- The coefficient of lateral force due 
to sideslip, Cy , was obtained from wind-t\innel data taken at M = O.I6 

for a 0.20-scale model of the test airplane, both complete and with tail 
removed. Changes in Cy due to increasing Mach nimiber were computed by 

applying equation (l) of^^reference 6 to the contribution of the vertical 
tail assuming that the tail-off value did not vary with Mach nimiber. 


Yawing moment due to sideslip .- The coefficient of yawing moment due 
to sideslip, Cn^^ was obtained from wind-tunnel data taken at M = O.I6 
and corrected for higher Mach numbers in a manner identical with that used 
for Cy • 

Rolling moment due to sideslip.- The coefficient of rolling moment 
due to sideslip, C^p, was determined from wlnd-t\innel data obtained at 
M = 0.16 at angles of attack of 0° and 8° for the complete 0.20-scale 
model, the model with tail removed, and the wing alone. 
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The wing contribution to 
follows : 

be I 




w 


SCl /m 


be I 


Pw 


SCt 


M=0 


was broken down into two parts as 


reference 6 


ac 




acL 


m=o.i6 

W.T. 


ac 




w 


ar 


acz 


Pw 


M 


ar /M=o 


reference 7 


ac 




w 


ar 


M=0.l6 

W.T. 


(Al) 


At Cl = 0^ the dihedral effect of the wing was assiomed to be due 
entirely to the geometric dihedral angle, reducing the first term of equa- 
tion (Al) to zero. The second term of equation (Al) was then assumed 
constant with Cl at a given Mach number, enabling the first term to be 
evaluated at lift coefficients greater than zero. The compressibility 
corrections indicated in equation (Al) were applied assuming that test 
results obtained at M = 0.l6 were essentially those at M = 0. 

The contribution of the vertical tail to was determined from 

the tail-on and tail-off data at angles of attack of 0^ and 8^ and cor- 
rected for higher Mach numbers using the method applied to the tail contri- 
bution to Cyp • 

The increment of C^p due to interference and the fuselage was 
obtained from wind-tunnel tests of the wing alone and the wing-fuselage 
combination at both 0^ and 8^ angles of attack and was assumed constant 
with Mach number. 


For the entire airplane, was determined in the following manner: 


CZp - (CZp)^ + (Czp)^ + (CZp)i+f 


(A2) 


Rotary Derivatives 


Rolling moment due to rolling .- The rolling-moment coefficient due to 
rolling velocity, was determined as a function of Mach number for the 
wing alone by application of figure 5 of reference 7* 
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The contribution of the horizontal tail was determined by the method 
applied to the wing. To express the values of for the horizontal 

tail in terms of wing area, span, and wing-tip helix angle, the values 
obtained from figure 5 of reference 7 were corrected by the following 
method : 

(Cn ) - 

^Clp^h - Sb2 

The contribution of the vertical tail was determined by the method 
presented in reference 3> using a value of sldewash parameter 

equal to 0.248, obtained from unpublished results of a theoretical investi- 
gation and tests conducted in the Langley stability tunnel. 

The contributions of the wing, horizontal tail, and vertical tail were 
added algebraically to obtain the estimated values of for the test 

airplane at different Mach nimibers and lift coefficients: 

Cip = (Cjp)^ + (CZp)^ + (Cip)y (a4) 

Yawing moment due to rolling .- The yawing-moment coefficient due to 
rolling velocity, Cnp, was determined as a function of lift coefficient 
according to equation (3I) and figure 25 of reference 8 for the wing alone. 
The variation thus obtained was corrected for compressibility effects by 
application of equation (3) of reference 6. 

The contribution of the vertical tail was determined by the method 
presented in reference 3 and added algebraically to the wing contribution 
to obtain Cnp for the entire airplane: 

Cnp = (Cnp)-vf + (Cnp)y (-^5) 

In the discussion it was noted that use of C^p and Cnp for the wing 
alone in the computations gave adequate agreement with experimental period 
and damping. 

Rolling moment due to yawing .- The variation of rolling-moment coef- 
ficient due to yawing velocity, C^^, with lift coefficient was determined 
from figure 26 and equation (37) of reference 8 for the wing alone. Cor- 
rections for compressibility effects were applied by means of equation (I5) 
of reference 6. 

The increment of Ci^ due to the vertical tail was fo\md by means of 
the following expression: 

2ivv. 

(^ir^v = ■ ^2“ a)(CYp)y 





p^h 


(A3) 




(a6) 
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For the airplane, is given as: 


C!r = (Clr)„ + (Clrt 


(A?) 


Yawing moment due to yawing .- The yawing-moment coefficient due to 
yawing velocity, was determined for the wing alone as a function of 

lift coefficient from figures 1 ^ and 27 and equation ( 4 l) of reference 8 . 

The contribution of the vertical tail was computed as a function of 
Mach number from the following relation: 


For each lift coefficient and Mach number considered, the estimated 
value of Cnp for the airplane is given as: 


The center of pressure of the vertical tail used to determine 
was obtained from figure 5 of reference 9 ^ using the aspect ratio, taper 
ratio, and sweep angle of the vertical tail, the root of which was assumed 
to lie on the fuselage reference axis. The center of pressure was assiomed 
to lie on the quarter-chord line. 

The lateral-force coefficients due to rolling and yawing velocities, 
CYp and were found to have little effect on the computed lateral 

motion of the test airplane. Therefore, those derivatives were assiomed 
to be equal to zero in this analysis. 



(A 8 ) 


Crir + (Cnr)v 


(A9) 
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TABLE I.- DIMENSIONS OF TEST AIRPLANE 


Wing 

Area 287-9 sq ft 

Span 37.12 ft 

Aspect ratio ..... ^-785 

Taper ratio O .513 

Dihedral . 30 

Sweepback of 0.25-chord line 35°l4* 

Root airfoil section (normal to 0.25-chord line) . . . NACA 0012-64 

(modified) 

Tip airfoil section (normal to 0.25-chord line) . . . NACA 0011-64 

(modified) 

Horizontal Tail 

Area 35 .O sq ft 

Span 12.75 ft 

Aspect ratio 4.65 

Taper ratio O. 45 O 

Sweepback of 0.25-chord line 3^°35’ 

Airfoil section (parallel to center line) NACA 0010-64 

Vertical Tail 

Area 39.75 sq ft 

Span 8.38 ft 

Aspect ratio 1-77 

Taper ratio 0-3% 

Sweepback of 0. 25 - chord line 350 


TABLE II,- PARAMETERS USED IN COMPUTING LATERAL PERIOD, DAMPING, AND JiEl 

|P| 


bp 

M 

Cl 

CLa 

°o 


Cnp 

CYg 

a 

^2p 

^Cip 

^Cnp 

Cnp 

CZr 

Cnr 




^1/2 

\\ 

|P| 

^Icpl 

|p| 

10,000 

0.35 

0.351 

0.0693 

-0.25 

-0.0836 

0.1060 

-0.680 

-0.3J^6 

-0.349c 

-0.0235 

-0.0308 

0.122 

-0.1770 

2.44 

2.43 

2.23 

2.30 

1.59 

1.58 


.211 

.0710 

-.40 

-.0684 

.1088 

-.686 

-.35*^ 

— 

-.0140 

— 

.097 

-.1795 

1.92 

— 

1.77 

— 

1.28 

— 

.55 

.142 

.0733 

-.40 

-.0590 

.1123 

-.696 

-.3620 

— 

-.0093 

— 

.086 

-.1840 

1.57 

— 

1.42 

— 

1.06 

— 

.65 

.101 

.0768 

-.40 

-.0565 

.1170 

-.708 

-.3700 

— 

-.0065 

— 

.080 

-.1887 

1.30 

— 

1.17 

— 

0.96 

— 

.75 

.076 

.0824 

-.40 

-.0565 

.1233 

-.724 

-.3800 

-.3860 

-.0049 

-.0026 

.077 

-.1953 

1.11 

1.11 

0.98 

0.98 

0.9c 

0.90 

35,000 

.55 

, hl 2 

.0733 

-.40 

-.0927 

.1123 

-.696 

-.3620 

-.3640 

-.0270 

-.0358 

.142 

-.1835 

2.56 

2.55 

3.85 

4.08 

2.01 

1.99 

.70 

.255 

.0791 

-.40 

-.0773 

.1199 

-.715 

-.3750 

— 

-.0160 

— 

.116 

-.1896 

2.02 

— 

3.06 

— 

1.61 

— 

.80 

.195 

.0869 

-.40 

-.0741 

.1273 

-.733 

-.3850 

— 

-.0120 

— 

.108 

-.1970 

1.73 

— 

2.63 

— 

1.45 

— 

.90 

.154 

.09*^3 

-.25 

-.0721 

.1366 

-.757 

-.3990 

— 

-.0092 

— 

.106 

-.2065 

1.49 

— 

2.21 

— 

1.31 

— 

1.00 

.125 

.0785 

-.50 

-.0768 

.1467 

-.782 

-.4l4o 

-.4190 

-.0068 

-.0047 

.104 

-.2170 

1.30 

1.30 

1.89 

1.87 

1.30 

1.23 


aValues coii5)uted for wing alone. 

bValues for wing plus vertical tail contribution computed by method of reference 3 * 


Constants: W = 12,500 Ix = 72^5 Cy^ = 0 

S = 287.9 Iz = 23,191 CYr = 0 

b = 37.12 e = 2.5 
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TABLE III.- AVERAGE VALUES OF THE MEASURED LATERAL OSCILLATORY CHARACTER- 
ISTICS OF TEST AIRPLANE AT VARIOUS MACH NUMBERS FOR ALTITUDES OF 10,000 
AND 35,000 FEET 


hp 

M 

p 

Ti /2 

ItI 

|P 1 

10,000 

o.Uo 

.50 

.60 

.70 

.79 

2.14 

1.70 

1.38 

1.15 
1.00 

2.05 

1.60 

1.35 

1.27 

1.21 

1.95 

1.70 

1.49 

1.30 

1.15 

35,000 

.50 

.60 

.70 

.80 

.90 

1.00 

l.ol^ 

2.55 

2.23 

1.95 

1.70 

1.50 

1.34 

1.27 

4.10” 

3.50 

2.92 

2.44 

1.85 

1.76 

2.00 

■ 2.48 
2.07 
1.91 
1.77 

1.64 

1.52 

1.48 
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Figiare 1 .- Photograph of the test airplane. 
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■e 3 .- Time histories of measured quantities during lateral oscil 
lation; Mach mmiber, 0.79; pressure altitude, 10,000 feet. 
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Figure U.- Time histories of measured quantities during lateral oscil- 
lation; Mach number, l.OU; pressure altitude, 35>000 feet. 
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Figure 5*- Variation with Mach number of experimental and computed lateral period and damping 

characteristics at each test altitude. 



ro 


NACA TN 3521 


NACA - Langley Field, Va. 





C3 


I 







Period, P, sec 


Figure 6 .- Comparison of experimental and computed lateral period-damping relationships with Armed 

Services specification dated 19 ^. 
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